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Abstract 

We study solutions to the free stochastic differential equation dXt = dSt — ^DV{Xt)dt, 
where y is a locally convex polynomial potential in m non-commuting variables. We show 
that for self-adjoint V, the law fiy of a stationary solution is the limit law of a random 
matrix model, in which an m-tuple of self-adjoint matrices are chosen according to the law 
exp(— A''Tr(V'(Ai, . . . , Am)))dAi ■ ■ ■ dA^- We show that ii V = Vp depends on complex param- 
eters Pi, . . . ,(3k, then the law nv is analytic in f3 at least for those /3 for which Vp is locally 
convex. In particular, this gives information on the region of convergence of the generating 
function for planar maps. 

We show that the solution dXt has nice convergence properties with respect to the operator 
norm. This allows us to derive several properties of C* and W* algebras generated by an m- 
tuple with law /xy. Among them is lack of projections, exactness, the Haagerup property, and 
embeddability into the ultrapower of the hyperfinite IIi factor. We show that the microstates 
free entropy xi'''v) is finite. 

A corollary of these results is the fact that the support of the law of any self-adjoint poly- 
nomial in Xi , . . . , Xn under the law fiy is connected, vastly generalizing the case of a single 
random matrix. 

1 Introduction 

There has been a great deal of interest in studying matrix integrals in physics since the work 
of 't Hooft who made the connection between the problem of enumerating maps and estimating 
integrals of the form 

Zj,{V) = y"e-^^(^(^--'^-))dXi ■■■dXm 

where dX denotes the Lebesgue measure on N x N Hermitian matrices and Tr the non-normalized 
trace Tr(^) = X^i^i^n- 1^ is a polynomial in m-indeterminates. Let us recall that a map of 
genus g is a graph which is embedded into a surface of genus g in such a way that the edges 
do not intersect and so that disecting the surface along the edges decomposes it into faces, each 
homeomorphic to a disk. We shall consider maps with colored edges and enumerate them when 
the degrees of the vertices, as well as the ditribution of color of the edges around each vertex, are 
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prescribed. The number of colors will be m, and the colors will be simply refered by the numbers 
{1, . . . ,m}. A vertex with colored half-cdgcs, an orientation and a distinguished half-edge, can be 
associated bijectivly with a non-commutative monomial q{X) = Xi^ ■ ■ ■ Xi^ as follows; the first (or 
distinguished) half-edge has color ii, second has color i2, etc., the last half-edge having color ip. 
Such a vertex, equipped with its colored half-edges, distinguished half-edge and orientation, will be 
called a star of type q. Wc will denote by Mg{{qi, ki)i<i<n) the number of maps with genus g and 
with ki stars of type qi for 1 < i < n, the maps being constructed by gluing pairwise the half-edges 
of the stars (the counting is done up to homomorphism of the surface and stars are labelled) . 
't Hooft showed that if 

^ m n 

v{x,, ...,x^) = ^(ft,,,),<,<„(Xi, ...,x^) = -J2xi + Y^Mx,,.. .,X^) 

i=l i=l 

where qi, i = 1, . . . ,n are monomials in m-non commutative indeterminates, then 

lU,(?ih<j<n g>Q fci,...,fc„eN" l<i<n * 

where the equality holds in the sense of formal series. Differentiating formally this equality, one 
also finds that if we consider the Gibbs measure 

then we have, for any monomial P, the formal expansion 



^^_^JP) := /lTV(P(X,,...,X„))rf^-^^^^^,^^^^jX^ 

= E]^ E n ^^M,((,„A;,)i<K„,(P,l)). 

g>0 ki,...,kn l<i<n * 

We refer the reader to the survey papers [13, 19] to see diverse uses of this fact in the physics 
literature. 

In [1, 12] for m = 1 and then in [20, 21, 22] for all m G N, these formal equalities were shown 
to hold in the sense of large N expansion when V satisfies some convexity hypothesis (or one adds 
a cutoff to make the integral finite) and the parameters (A)i<i<n are chosen to be small enough. 
In [20] , one of the key steps towards this analysis is to notice that converges towards 

W'9i>l<i<n 

a limit, denoted later ^ > which satisfies the so-called Schwinger-Dyson equation 

T®T{diP)=T{DiVP) (2) 

for all polynomials P and all i G {l,...,m}, and with V = W^(ft.(jj)i<i<„- Here, di and Di are 
respectively the non-commutative derivative and the cyclic derivative with respect to the variable 
Xi (sec the next section for a definition). One then shows that for sufficiently small parameters 
this equation has a unique solution, which is exactly the generating function for planar maps: 

^^(/3.,,),,.,„(^)=^(ft,.0i<.<n(^)^= E n ^^Mo((gi,fcOl<Kn,(P,l)). 

ki,...,kn l<i<n 



2 



It is natural to wonder how to use these expansions to study the numbers Mg{{qi, fci)i<i<„) and 
in particular their asymptotics as A;i , . . . , fc„ go to infinity. The answer to this question is still open 
in such a general context. It is, however, quite well understood in the case m = 1. Let us highlight 
this point in the case of quadr angulations, corresponding to the potential V{x) = Px^, /3 G M, even 
though the enumeration of quadrangulations was achieved by direct combinatorial arguments by 
Tutte [23] long ago. In this case, the Schwinger- Dyson equation (2), taken at P{x) = {z — x)~^ 
shows that the Cauchy transform 

J Z — X P'"^ J z — X 

satisfies an algebraic equation of degree two 

/3 3 
^ ~ ^ dTw ^4 (X) - 1. 
X Z ' 

We can solve this equation in terms of P{z) which is a polynomial of degree two with two unknown 
coefficients; we then find that G is given by a polynomial plus the square root of a polynomial of 
degree six, that we denote Q, with two unknown coefficients. Until this point, all the arguments 
follow the induction relations already found by Tutte. However, the difference now is that we know 
that Mis^x"^ = Tw^ 4 is a probability measure on M. Assume we can argue that for sufficiently small 
P, the support of tw^ ^4 is connected. Then, we see that, because this means that G is analytic 
outside an interval, the polynomial Q must have two double roots. This actually determines Q, 
and thus G, uniquely. Since G is the generating function for planar maps, we are done. Hence, we 
see in this context that the a priori information that -^(flj.gj)i<;<„ is the Cauchy transform of 
a measure on the real line (which is not clear from its definition as a generating function of maps), 
and with a connected support, is enough to conclude. 

The goal of this article is to push forward the analysis of the limiting tracial state Tw,n , 
in the multi-matrix context. We prove in particular that when I^(/3j,ijj)i<j<„ satisfies a certain local 
convexity property (see (3)), the support of the limiting spectral measure of the random matrices 
with law u,^ is connected. In fact, the same is true for an arbitrary non-commutative 

polynomial in the random matrices. Note that this information is enough to solve the enumeration 
problem when m = 1 as we have seen above for quadrangulations (though connectivity of support 
can in this case be proved by other techniques, see e.g [10]). 

The tracial states r = ry under consideration will be solution to the Schwinger-Dyson equation 
(2) for some general potential V. Non-commutative laws arising as limits of laws of random matrix 
models given by (1) have also naturally appeared in free probability theory. There, the fact that 
they satisfy the Schwinger-Dyson type equation is restated as the fact that the free conjugate 
variables of the law are equal to the cyclic gradient of a polynomial potential, see also [8, 3]. 

In the multi-matrix setting, uniqueness of the solution to the Schwinger-Dyson equation is 
unclear in general. It was proved in [20] that, when V = ^^(^i,qi)i<j<„) there exists a unique 
solution such that |T(Xjj •••Xj^)| < for all choices of ij G {!,••• ,m} and all k = 1,2,..., 
provided the PiS are sufficiently small. In this paper we define a notion of locally strictly convex 
potential, which generalizes to non-commutative variables the standard notion of local convexity 
for functions on the real line. One of the central result of this paper will be the uniqueness of the 
solution to Schwinger-Dyson equation for locally strictly convex potentials V, when the domain of 
strict convexity is large enough. 
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We also show that if F = Vjg = ^ Piqi with some monomials qi and (3 = (A)i<i<n) P TVp{P) 
is analytic in the whole region of local convexity of Vg, for any monomial P. Because 

ki,...,kn l<i<n * 

this result shows that there is no breaking of analyticity of (/^i)i<2<n 

domain where l^(gj,/3j)i<j<„ stays locally strictly convex, and thus provides valuable information on 
the asymptotics of the numbers Mo{{qi, ki)i<i<n, {P, 1)). 

Of particular interest are the C*-algcbras {Av,tv) and the VF*-algebras {My,tv) generated 
by operators having a law satisfying the Schwinger-Dyson equation with a fixed locally convex 
potential V. We derive several properties of such algebras, showing that they are projectionless, 
exact [28] and possess the compact approximation property of Haagerup [16], and can be embedded 
into the ultrapower of the hyperfinite IIi factor. We also show that the algebras My are factors 
and that the generating operators have a law with finite free entropy. This has as consequences a 
number of properties of the algebras My (including primeness and lack of Cartan subalgebras) , see 
[26] and [15]. All of these properties are similar to (and are often derived from) the corresponding 
properties of free group factors (which are von Neumann algebras associated to a quadratic potential 
V{Xi, . . . , Xn) = ^ Y2i ^i)- This adds evidence towards a positive answer to Voiculescu's question 
of whether My is isomorphic to a free group factor for a fairly arbitrary potential V. 

While somewhat technical, it should be noted that studying properties of Ay and My is of 
substantial interest. Indeed, connectivity of support of limit distributions of random matrix models 
is directly related to the lack of projections in the C*-algebra Ay (which in turn is derived from 
the famous result of Pimsner and Voiculescu [24] , essentially dealing with the quadratic V) . 

Finally, in the remaining sections of the paper we show that the random matrices following the 
law (1) give a very good approximation to the non-commutative law Ty. If V is locally convex, we 
show that the lim sup in the definition of the microstates free entropy of Ty [25] can be replaced by 
a limit. In the case that V is (globally) convex, we show that operator norms of arbitrary polyno- 
mials in such random matrices almost surely approximate the operator norms of such polynomials 
computed in the C*-algebra Ay. This extends the results of [14], which correspond to the case of 
quadratic V (our proof, though, relies on their result). 

The main technical tools used in the present paper involve the study of free Langevin-type 
diffusion and its convergence to a stationary measure which corresponds to the limit law Ty of 
random matrices following the measure (1). We extend some of the results of [5] to the setting 
of locally convex potentials (sec below). This way, we are able to show that operators having 
a specific limit law can be approximated in the operator norm by continuous functions of free 
Brownian motion. This enables us to carry over a number of properties of the algebra generated 
by free Brownian motion to the algebras Ay and My. 

1.1 Definitions, notations and statment of the results. 

Let us now state more precisely our setup and results. 

We let C(Xi, . . . , Xm) be the set of polynomials in m non-commutative variables (Xi, . . . , X^)- 
We shall not assume in general that {Xi, ... , X^) are self-adjoint but let {X^, . . . , X^) be their 
adjoints for some involution *. We denote C{Xi, . . . ,Xjn,X*, . . . ,X^) the set of polynomials in 
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the non-commutative variables {Xi, . . . , X^, X^, . . . ,X^). This set is endowed with the linear 
involution so that 

\ h ik ' ik ik-i H 

where wc denoted in short X} = Xi and X^^ = X* and the (ei, . . . ,ek) belong to {—1, 1}^. We 
shall denote below for two sets of non-commutative variables {Xi, . . . and (Yi, . . . and 

an involution * 

1=1 

II • lloo will denote an operator norm such that the completion of {C{Xi, . . . ,Xjn,Xi, . . . ,X^), *) 
for this norm is a C*-algebra. For a m-dimensional vector X = {Xi, . . . ,Xjn) we denote in short 

^lloo = Ill&Xi<j<„j 1 1 Xj lloo- 

We let V G C{Xi, . . . , X^) be a polynomial in m non-commutative variables. We will say 
that V is (c, M) convex iff for any m-tuples of non-commutative variables X = (Xi, . . . and 
Y = {Yi, . . . ,Yjn) in some C*-algebra {A, \\ ■ ||oo) satisfying ||Xj||oo, ||^||oo < M, z = 1, . . . ,m we 
have 

[DV{X) - DV{Y)].{X -Y)> c{X - Y).{X - Y) (3) 

where the inequality is understood in the sense of operators {X >Y iS X — Y is self adjoint and 
has non negative spectrum). D = {Di, . . . ,Dm) denotes the cyclic gradient which is linear and 
given, for any monomial P, by 

DiP = J2 ^2^1. 

P=PlX,P2 

Later, we shall also need the non-commutative gradient d which is given, for any monomial P, any 
i G {1, . . . ,m} by 

diP = ^ Pi ® P2. 

P=PlXiP2 

We occasionally shall consider polynomials in C(Xi, . . . , X„i, X^, . . . , X*^); in that case we extend 
di and Di by setting diXj = lt=jl (X" 1 and diXj = ^ whereas we have also the derivative 
and Di^^ with respect to X* which satisfy di^^Xj = (8) but di^^X* = 1 (g) 1, and extending by 
linearity and Leibnitz rule. 

In the case that Xi, . . . ,X„ arc self-adjoint, we shall make the convention that diXj = lj=jl®l, 
while di^^^Xj = for all j (in other words, we shall continue to think of all quantities as functions 
of Xi , . . . , Xn alone) . 

Note that by taking X = Y+eZ and letting e going to zero we find that for any bounded operator 
Z and any operator Y with norm strictly smaller than M, the condition that V is (c, M)-convex 
implies that 

in rn 

YI {diDjViY)^Zi X Z* + Zj X idiDjViY)^Zi)*) > 2cZ.Z. 

1=1 j=i 

Wc shall say that V G C{Xi,...,Xm) is self-adjoint iff for any self-adjoint variables X = 
(X, • • • , Xm), [ViXi, Xm)]* = F(Xi, . . . , X^). V is self-adjoint (c, M)-convex if V is self- 
adjoint and the above holds once restricted to self-adjoint variables, i.e for any m-tuples of self- 
adjoint variables X = (Xi, . . . ,X^) and Y = (Yi, . . . ,1^) living in some C*-algebra {A, \\ ■ ||oo) 
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which are bounded in norm by M, we have 



[DV{X) - DV{Y)].{X -Y)> c{X - Y).{X - Y) 

In this case X.Y = {X,Y} := ^ + YiXi) is simply the anti-commutator of X and Y . 

If we speciahze this assumption to matrices and consider A to be the algebra oi N x N matrices 
with complex entries equipped with the usual involution {A*)ij = Aji and the spectral norm || ■ ||oo, 
we find that if V is self-adjoint (c, M)-convex, {Xij)i<j Tr[V{X)] is strictly convex on the set of 
entries where X is Hermitian and with spectral radius bounded by M since 

TtV{X) -TtViY) = [ TtiDViaX + {l-a)Y).{X -Y))da 
Jo 

= Tr{DV{Y).{X -Y))+ f Tr{[DV{aX + (1 - a)Y) - DV{X)].{X - Y))da 

Jo 

> Tt{DV{Y).{X - Y)) + ^Tt{{X - Yf). 

Taking Y = {X + Z)/2 and X to he X or Z and summing the resulting inequalities gives 

(X + Z\ 
— ^ j > cTr((X - Zf) 

and hence the Hessian of {Xij)i<j TrV{X) is bounded below by cl, at least on matrices X with 
norm bounded by M. 

This kind of hypothesis was shown to be very useful in [20]. The interest in relaxing the 
hypothesis of convexity to hold in a bounded domain is related with matrix models where V = 
^X.X + W with W = J^i^i iSiQi for some monomials and complex parameters {Pi)i<i<n G C". It is 
clear now that for all M finite we can choose the /3j's sufficiently small so that V is (1/2, M)-convex 
(whereas it would not work with no bounds). Indeed, in that case 

{DV{X) - DV{Y)).{X -Y) = {X- Y).{X -Y) + {DW{X) - DW{Y)).{X - Y) 

But when the norms of X and Y are bounded by M, 

\{DW{X) - DW{Y)).{X -Y)\< C{M) max |A|(X - Y).{X - Y) 

i 

with C{M) a constant which only depends on M and the qi. Hence, we can now choose t small 
enough so that C{M) maxj \ti\ < 1/2 and so V is then (1/2, M)-convex. This is analogous to what 
was done in [20] in case of non-convex interaction; it was shown that then if one adds a cut-off, the 
large N expansion is still valid provided the parameters in W are small enough. 

Hereafter we assume that V is (c, M)-convcx. Wc let (A, *,([>) be a non- commutative probability 
space generated by a free Brownian motion S (wc refer to [5, 4] for an introduction to free Brownian 
motion and its related free Ito calculus). We shall denote by || • ||oo the operator norm in (A, </>). 

We prove (see Lemma 2.1 and Theorem 2.2): 

Theorem 1.1. Let V be a (c, M) -convex polynomial in Xi, . . . , X^- 

Then there exist Mq = Mo{c,\\DV{0).DV{0)\\oo), Bo = Bo{c,\\DV{0).DV{0)\\oo), and b = 
b{c, \\DV{0).DV{0)\\oo, M) > Bq finite constants, so that whenever M > Mq and Z is an m-tuple 
with \\Z\\ < b, 
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i. There exists a unique solution Xf to 



dXf = dSt - ^DV{Xf)dt, t e [0, +00), (4) 



with the initial data X^ = Z. Moreover, in this case, 



\\x^\\oo<M, Vie [0, +00), 

limsup lloo < -Bo, 

t—*oo 

Xf G C*{Z, Sq-.qe [0, t]), Vt G [0, +00). 

a. \\Xf - X^Woo -^0 as t-^ 00. 

Hi. The law of {{X^)*,X^) converges to a stationary law fiy £ C(Xi, . . . , Xm, X*, . . . , X^y as 
t goes to infinity, fiy is the non- commutative law of m variables uniformly bounded by Bq. If 
V is self-adjoint, nv is the law of m self-adjointvariables. 

iv. The restriction Ty = /Uy|c(Xi x,n) of fiy to C{Xi, . . . , Xm) satisfies the Schwinger- Dyson 
equation which states that for all P G C(Xi, . . . , X„i) and a// i G {1, . . . , m} 

Tv®rv{diP)=Tv{DiVP). (5) 

Moreover, if v is the law ofm variables whose uniform norm is bounded by b and i'\c{Xi,...,Xm.) 
satisfies (5), then iy\c(Xi,...,Xm) = '^v- 

The Schwinger-Dyson equation (5) is exactly the same as the one which characterized the 
enumeration of maps (2) from which we deduce that ■A^(^j,qj)i<j<„ = ri x^+J2P-q- least for 
sufficiently small polynomials V. This allows to give some information on the domain of analyticity 

of (3= (A)i<.<n ^ -^(ft,g,)i<i<„(^) = Mp{P) (see section 3). 

Theorem 1.2. Let V = = Y17=i f^ili ^ polynomial, where fi = {Pi)i<i<n are (complex) 
parameters and {qi)i<i<n are monomials. For c,M positive real numbers and M > Mq with Mq as 
in Theorem 1.1, let T{c,M) C C" be the interior of the subset of parameters /? = {Pi)i<i<n such 
that Vfs is {c, M) -convex. 

Then , for any P G C{Xi, . . . , Xm) , (3 G T{c,M) tv^{P) is analytic. In particular, (3 — >■ 
JV\.j)[P) extends analytically to the interior of the set of Pi's where ^Yl^Li^f + {c,M)- 
convex for M > Mq (c) . 

The laws fiy are interesting in their own. In the free probability language we have proved the 
following. 

Theorem 1.3. Let V be a {c, M) -convex potential with M > Mq the constant of Theorem 1.1. IfV 
is self-adjoint, there exists a non- commutative law in C{Xi, . . . , Xm, X^, . . . , X^)' with conjuguate 
variable {DiV)i<i<,m- There exists at most one such law satisfying the additional constraint to 
be the law of variables bounded by b. For non self-adjoint potential, there a unique law nv £ 
C{Xi, . . . , Xm,Xl, ... , Xm)' which satisfies for all P G C(Xi, . . . , Xm, XI,... , Xm) and all i G 
{1, . . . ,m} 
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m m 

^ywv ® {{di + 5i,,)(A + Di,,)P) = ^ /xy {DiVDiP + {DiVyDi,,P) (6) 

1=1 i=l 
with the non- commutative deriavtives with respect to X* . There exists at most one such 

law satisfying the additional constraint to be the law of variables bounded by b. 

This statement is deduced from (2.2) (2). Moreover, we let Z be an m-tuple of operators with 
law We shall prove that the C*-algebra and the von Neumann algebra generated by Z have 
many properties in common with the one generated by a semi-circular system. 

Theorem 1.4. Assume that V is {c, M)-convex with M > Mq the constant of Theorem 1.1. The 
C* -algebra generated by Z is exact, projectionless and its associated von Neumann algebra has the 
Haagerup approximation property and admits and embedding into the ultrapower of the hyperfinite 
III factor. 

In particular wc have 

Corollary 1.5. Assume that V is {c, M)-convex with M > Mq the constant of Theorem 1.1. Let 
Z be any m-tuple of b-bounded variables with the unique law jiy satisfying (2). (a) The algebra 
C*{Z) has no non-trivial projections, (b) The spectrum of any non-commutative *-polynomial P 
in the m-tuple Z is connected (in the case that P{Z) is normal, this means that the support of its 
spectral measure is connected), (c) If P is any polynomial in Z whose value is self-adjoint, then 
the probability measure given by the law of P{Z) has connected support. 

2 Existence of free diffusions and convergence to their stationary 
measure 

We shall show that if M is chosen large enough (depending on c and llL'y(O). 1)^(0) ||oo), we can 
build a bounded solution to some free stochastic differential equation with drift DV provided that 
V is (c, M)-convex. This generalizes Langevin dynamics to the context of operators. 

Lemma 2.1. Let V be a (c, M)-convex polynomial in Xi,. . . , Xm. Then there exist finite constants 

Mo = Mo(c, II W(0).W(0)||oo), Bq = Bo{c, || W(0).W(0)||oo) 
b = b{c,\\DV{0).DV{0)\\oo,M)>Bo 

so that if M > Mq, and Z is any m-tuple with \\Z\\ < b, there exists a unique solution Xf to 

dXt = dSt - ^DV{Xt)dt, t e [0, +oo). (7) 

with the initial data Xq = Z. Moreover, in this case, 

||Xt||<M, VtG[0,+oo), 
limsup \\Xt\\oo < Bq, 

XteC*{Z,Sq:qe[0,+oo)), Vt€[0,+oo). 

IfV is self-adjoint (c, M) -convex, the above results hold under the additional assumption that 
Xq = Z is self-adjoint. In this case, Xt remains self-adjoint for all t > 0. 
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Proof. We remind the reader (cf. [5]) that if DV is uniformly Lipschitz for the uniform norm, 
the existence and uniqueness to (7) is clear by the following Picard argument. For the existence 
we consider the sequence X", n = 0, 1, 2, . . . constructed recursively as follows. Set = Z for all 
t, and having defined X", let X""*"^ be given by the equation 

dX^+^ = dSt - ^DV{X^)dt 

with the initial condition Xq^^ = Z. 

Subtracting the equations for X^~^^ and X" from each other, we get for all t > 0, 

\\X^+' - X-\\oo <\f^ \\DV{X-) - DV{Xr')\\oods < \\\DV\\c \\X- - X^'' Woods 

where ||-Dy||£ denotes the Lipschitz norm of DV 

\\DV\\c = inf{C : \\DV(Y) - £>F(y')lloc < C\\Y - Y'\\^}. 
Iterating, we deduce that 

-sup||5„ + uW(Z)|U< _(2Vt + -||W(Z)|U) 

which proves norm convergence of X" (note that Xq = Z for all n). Moreover, this limit satisfies 
(7). Wc also see that G C*{Z, Sq : q > 0) and therefore the limit Xt G C*{Z, Sq : q > 0) as well. 
The proof of uniqueness follows the same lines since any two solutions Xt,Yt satisfy 

||Xt-Yt||oo< \\Xs-Ys\\oods 

which proves that X = y by Gronwall's argument. 

If y is a self-adjoint polynomial, then D^V is also self-adjoint for all i G {1, . . . ,m}. Therefore, 
since St is self-adjoint for all t > 0, we deduce by induction that X" is self-adjoint for all n > 
and all t > and so its limit Xt is also self-adjoint. Hence the solution Xt of the free stochastic 
equation is self-adjoint for all times. 

We now return to the case of a (c, M) convex V. By Lemma 3.2 in [5], we can construct a new 
function 

fR{x) = ^DV{x)h{j2 mu 

i 

so that fji is uniformly Lipschitz and //j(X) = DV{X) if ^ ||Xj||oo < R- The Picard argument 
above implies existence and uniqueness of a solution Xf' to 

dXf = dSt + fR{X^)dt 

for all times t. Clearly, if we show that this solution satisfies J2j ll(^/)t||oo < R, it will also be a 
solution to the original equation (7) involving DV. Furthermore, if we start with some initial data 

Xq with ^ ||(Xj)o||oo < R, solutions to (7) always exist for small time (at least up until the time 
that the operator norm of the solution exceeds R) . Thus we may consider a solution up to the time 
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that its norm reaches some fixed constant M (with the intent to show that this time is infinite). 
By free Ito calculus 



dXt.Xt = 2Xt.dSt- DV{Xt).Xtdt + 2dt 

= 2Xt.dSt - {DV{Xt) - DV{0)).Xtdt - DV{0).Xtdt + 2dt 

Therefore, for all s > 0, 

e'^'Xs.X, = Xo.Xo-2 r e''DV{0).Xtdt- r e'''[{DV{Xt) - DV{0)).Xt - cXfXtjdt 

Jo Jo 

+2 / e'^dt + 2 I e^KdSt 
Jo Jo 



< Xq.Xq 



[ e^L>F(0).. 
Jo 



Xtdt + 2c-^e'^ + 2 / e'^.dSt 



f 

Jo 



„ct 



where the last inequality holds in the sense of operator and we relied on our hypothesis of (c, M) 
convexity. Since also e'^^Xg-Xg is a non negative operator, we deduce that 



-^5 • — ||-^S '-^S II OO 



< e-^'Ao + 2c~^ + 2 



r 

Jo 



e<^-'^Xt.dS-, 



1 /"S 



+ ||W(0).W(0)||^ / e^^^-'^Afdt 



By Theorem 3.2.1 of [4], we know that the free analog of the Burkholder-Davis inequality for 
integrals with respect to free Brownian motion holds for the norm even with p = oo. More 
precisely, the following estimate holds: 



e^^^-'^Xt.dS, 



<2V2{ I e^'^^-'^Atdt] <^^yj e<^-'^Atdt 



Moreover, by the Cauchy-Schwarz inequality, we obtain the bound 



a 



\<t-^)A}dt 



' lo ■ <\f^ e<^-'^Atdt. 

1 

Hence, we get the inequality (with C = ^{\\DV{0).DV{0)\\io + 4v^)2 and since < b'^rn) 

Al < 4e-^'m%^ + 2'^c-^ + C r e<'-'^Atdt 

Jo 

< 4e-'^^mV + 2V2 + - f e<*-') {B'^ Al + B)dt 

< 4e-'=^m26^ + 2^0'^ + CB{2c)-^ + Cfi"^ /' e<^-''>Ajdt 

Jo 

where B is any positive constant (we just used that for all x, 2x < B + B~'^x'^). 

We now use Gronwall's lemma (or simply iterate in the above inequality) to deduce, with 
C" = 2^0-2 + C5(c)-S that 

Ay < [Am%'^ + C'e"'] + % f e^^'-^'^iWb^ + C'e^''\du 

B Jo 
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We now take B large enough so that c > ^ to conclude that 

A'i<8mWi-''> + C'— 



c- ^ 

^ B 



If we now choose Bq = ^ + CB{c) ^)-fu; for B = 2C/c, we have shown that: 

• stays bounded by 8m^6^ + B^. So if 8m^6^ < — Bq, Ag always stays bounded by M. 

• As s goes to infinity, limsupg^o^ As is bounded by Bq. 



• We can choose 6 > i?o as long as M > Sqa/I + Sui'^Bq := Mq. 

• XteC*{Z,Sg:q>0). 

This concludes the proof. □ 
We now consider the solutions of Lemma 2.1 starting with different initial data. 

Theorem 2.2. Let Mq, Bq and b be as in Lemma 2.1, and assume that M > Mq, and that Z is 
an m-tuple of operators with ||2'||oo < b. Consider the unique solutions Xf , to the free SDE 

dXt = dSt - ^DV{Xt)dt (8) 

with initial conditions Xq = Z, Xq = 0. Then 

1. \\Xf - X^W^ ^0 ast^oo. 

2. The law of {{Xf )*,Xf ) converges to a stationary law jiy £ . . . , X^, . . . , Xj^)' 
which satisfies for all P G C{Xi , Xm, Xf,..., X^) and all i e {1, . . . , m} 

m m 

Y,l^y®l^y + 9i,*){Di + A,*)-P) = ^l^v {DiVDiP + {DiVfDi^.P) . (9) 
i=i i=i 

Moreover, for any A; G N, any z G {1, . . . , m}, 



[{x,xtf) < BlK (10) 



Any law v G C{Xi, . . . , X^, X\., . . . , X^)' of variables bounded in operator norm by b which 
satisfies (9) equals /j, on C{Xi, . . . , X^, X^, . . . , X^) . 

3. The restriction Ty = A*i^|c(Xi Xm) of nv to C(Xi,...,X^) satisfies for all polynomials 

PeC{Xi,...,Xm) 

m m 

^ ry ® MdiDiP) = J2 MDiVDiP). (11) 

i=l i=l 

Any law v G C{Xi, . . . , X^n, X^, . . . , X*y of variables bounded in operator norm by b whose 
retriction iy\c{Xi,...,Xm) satisfies (11) is such that iy\c{Xi,...,Xm) = '^v- 
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Remark 2.3. (a) Note that if Z is an m-tuple of variables, norm bounded by b having the stationary 
law fi, then is a stationary process (with law given at all times by /i) and Xf — X^ converges 
to zero in operator norm as t goes to infinity by the first point of the theorem. 

(b) Recall that a law fi oi m non-commutative variables has conjugate variables (^i)i<i<m if 
and only if for any P e C(Xi, . . . , X^, X* , . . . , X;;) , 

Taking the adjoint, we find that we must also have 

ti®ti{di,.P)=lx{CiP)- 

Taking P = DiQ in the first equality, P = Di ^^Q in the second and summing the resulting equalities 
yields 

H (g) ii{{diDi + = ^l{i^DiQ + C*Di,*Q) 

which differs from (9) by the terms fi ® ^{diDi^^Q + di^^D-iQ). Hence, (9) is not equivalent with 
the fact that has conjuguate the variable {DiV)i<i<m in the case that V is not self-adjoint. In 
the self-adjoint case, because of our convention, Q depends only on Xi, . . . ,X„ and so the terms 
involving Dj^* and Sj,* are equal to zero. In that case, (9) is compatible with the condition that 
the conjugate variables are equal to the cyclic gradient of V . 

Proof. Consider two solutions Xj , Xf with initial data Y and Z, respectively, and assume 
that ||y|| oo, ||-^||oo ^ n^- Then 

d(Xf - XY) = -^[W(Xf ) - DV{Xr)]dt 

Since by Lemma 2.1, the operator norms of Xf and X^ stay bounded by M for all t and V is 
(c, M) convex, we find that 

d(Xf -Xr).(Xf -Xf) < -[DV{Xf)-DV{X}')UXf-XY)dt 

< -c{xf - xr).{xf - xr)dt 

where the inequality again holds in the sense of operators. This implies that 

||(Xf - xYUXf - Xr)||oo < e-^\\Z - Y\\l (12) 

and so 

hm ||(Xf-Xr).(Xf-Xr)||oo = 0. 

t— >oo 

In particular, we can take Y = and we have proved that all diffusion solutions starting from 
different initial data norm-bounded by b will asymptotically be the same as X^. 

As t gets large, X^ is bounded hy Bq < b according to Lemma 2.1, and so we can choose Y = Xf 
(with X constructed as X but with some free Brownian motion S on [0, s] and the increments of 
5 on [s, s -I- ) to deduce 

||(X0 - X^VJ.(X0 - Xo+Jlloo < e-'^'M. (13) 

Since X^^^ has the same law as X^^s we conclude that the joint law of X^, (Xf)* converges as t 
goes to infinity. We denote this limit by /i. Clearly, is a stationary law for the diffusion (since 
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-^oo+s ^'^d X^^^ have the same law). We now consider X^° to be the solution of the free SDE 
starting from X so that {X,X*) has law fi. Noting that we have 

dX^ = dSt - ^DV{X^)dt d{X^y = dSt - ^{DV{X^)ydt 
and applying free Ito's calculus [4], we have that for any P € C(^i, . . . , Xjji, X*, . . . , X^) , 

o = dtcp{p{xr,ixn*)) 

= ^^4>{^ + di){Di,* + Di)P{X^, {Xn*)) 

i=l 

- ^cf>{DViXn.DP{Xn + {DV{Xnr.D.P{X^, (^D*))) 

so that for all *-polynomials P G C{Xi, . . . , X^, X^, . . . , X^), jiy must satisfy (9). tv, the restric- 
tion of ixv to C(Xi, . . . must then satisfy (11). 

The uniqueness of solutions to this equation is simply due to the fact that if we run the process 
from Z bounded uniformly by b with law v satisfying (9), the law of X^ must be stationary, 
but also converging to fiy by the previous argument. Hence it must be equal to jiy. The same 
argument applies in the case that we only care about the restriction to C(Xi, . . . of some 

law V of variables Z bounded by 6, since if it satisfies (11), the process X^ will be such that 
dt(t>{P{Xf)) = for all P G C{Xi, X^) and so the law of Xf restricted to C(Xi, . . . , X„) will 
be stationary (here we use the fact that DV{X) depends only on X and not on its adjoint when 
speaking of the evolution of the law of Xf restricted to C(Xi, . . . , X^)). Since it also converges to 
Ty, we obtain the desired equality. □ 

We shall see in section 3, Corollary 3.2, that in fact /x not only satisfies (11) but actually for all 
zG{l,...,m}, allPGC(Xi,...,X^), 

li®li{diP) = ii{DiVP). 

In other words, at least when /x is the law of self-adjoint operators, the conjugate variables of /x are 
in the cyclic gradient space. 

3 Analyticity of the solution to Schwinger-Dyson equation and 
discussion around phase transition 

We show in this section that on the domain where V stays (c, M)-convex for some c > and 
M > Mo as in Lemma 2.1, the law ry will depend analytically on the parameters of V. 

Lemma 3.1. Let V = Vp = PiQi be a polynomial, where {3 = (/3i)i<i<n o-i"e (complex) 

parameters and {qi)i<i<n 0,1"^ monomials. For c,M positive real numbers and M > Mq with Mq 
as in Lemma 2.1 and Theorem 2.2, let T{c,M) C C" be the interior of the subset of parameters 
(3 = {Pi)i<i<n for which V is {c, M)-convex. Let iip = i^Vp be the unique stationary measure of 
Theorem 2.2 and rg the law of {Xi, • • • , X^) under fip. 

Then for any polynomial P G C{Xi, . . . , Xm), the map (3 G r(c, M) — ^ Tp{P) is analytic. 
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Note that T{c,M) is non empty as soon as the set of monomials {qi,l < i < n) contains 



{Xl,l<i<m). Indeed, if we set V{X) = YT=\ Pi^f + ELm+i we always have that for 



X, Y uniformly bounded by M, 

/ n n \ 

V PiDqiiX)- V PiDqi{Y)].{X-Y)<C{M) max mX-Y).{X-Y) 



\i=m+l i=m+l 



with a universal constant C(M) which only depends on M and the {qi)i<i<n- Hence, if /Jj > for 
ie {!,■■■ ,m}, 

{DV{X)-DV(Y)).{X -Y)>[uim I3i-C{M) max \Pi\]{X -Y).{X -Y). 

i<i<m m+l<i<n 

Therefore any set of parameters (A)i<i<n such that 

min Pi - C{M) max > c 

l<i<m m+l<i<n 

will be such that V is (c, M)-convex. 
Proof. 

We denote (xf )t>o the solution of (8) with potential V = Vp and starting from the null operator. 
We shall show that (5 —>■ X^ expands as a sum of uniformly bounded operators. More precisely, we 

(k-i k ) 

fix /? in the interior of T(c, M) and find a family "' " , fc^ G N, 1 < i < n of operator-valued 

processes such that for rj G C"^, \(3 — r]\ := maxi<j<„ — r]i\ small enough, 

n 

X^ = X,^+ U(^^-P^)''Xi'^-'-^ (14) 

ki fc,i6N'i j=l 

Moreover, X^ are operator-valued processes such that there exists a constant C which only 

depends on c, M and the degree of V so that 

sup ||xf ''•••'''"^lloo < C^'''- (15) 

Finally the distribution of 

^^^(fei,...,fc„)^^ ^ converges (in the sense of finite marginals, i.e., 
on polynomials involving only a finite number of the ( - ) ) towards the law of 

( X^^''"'''"^] as t goes to infinity. 

V / fei,...,fc„GN" 

Let us conclude the proof of the lemma assuming (14) and (15). (14) and (15) entail that for 
all polynomial functions P G C{Xi, . . . , Xm), for all t > 0, 

P ^ </>(P(xf )) 

is analytic in the interior of T(c, M) since it implies that for /3 G T(c, M), 

MP{x?)) = ct> Ip{x^ + fl{vi-Pi)'"xi'^'-''-^) 
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for all ri in the domain B{C,P) = {\r] — P\ < l/C*} which does not depend on the time parameter 
t > 0. Note also that t{P{X^)) is uniformly bounded independently of i G M+ since C does not 
depend on t. We know that (X^) converges as t goes to infinity towards 

n 

(note here that convergence of the X^'^'^''"''^'") in the sense of finite marginals is sufficient since 
^fei k„-Y,ki>K nr=i('?i — Pi)^^X^'"''^" goes uniformly to zero as K goes to infinity). But then the 
limit has to depend analytically on /3 (as a limit of uniformly bounded functions which are analytic 
on a fixed domain). This proves the claim. 

We now prove (14) and (15). We first check that (3 X^ is of class C°°, then that it is 
in fact an entire function and bound uniformly its radius of convergence. Finally, we prove that 
(^f'"'"''"'')fci,...,fc„eN« converges. 

Step 1: Pe T{c, M) is of class for aUt>0. 

Let us study the first order differentiability, and first check that /? — ^ X^ is continuous. In fact, 
if lp(i) = for i / p and Ipip) = 1, we write 

xf+^^^'-Xf = IjDVfs+^^^ixt^'n - DV^{X^)]ds 

= -\ft 5Z)F^+ei,((l - oc)Xl + aXf+^'-)tt[xf+^'- - Xl]dads 
^ Jo Jo 

~\L [DVi5+ei,-DVp\{X^).dt 

We find that if /3 + elp and (3 both belong to r(c, M) so that Xl^ and X^+'^^p stay uniformly 
bounded by M, that 

Jo 

where C only depends on (c, M). So Gronwall's lemma shows that for any time t >0, there exists 
a finite C{t) {C{t) is uniformly bounded on compacts) so that 

||xf+^'''-Xf||oo<C(i)e. (16) 

This suggests that /3 xf is in fact differentiable. To prove this point, let us introduce the 
candidate for the corresponding gradient; we define (V/jXf '')i<j<OT to be the m-tuple of operators 
valued processes solution of 

m 

dV^xf = - ^ d,D,V{X^)^V^pX^'Ut + {dFpD,Vp){X^)dt V^Xo^'^ ^ 

where dp is the standard gradient with respect to the parameters (3 (so d-'^DiV^ = Diqj for j G 
{1, . . . , n}). Here p is any integer in {1, . . . , n}. There is a unique solution to this equation (since 
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it is a linear differential equation as is given). By the same type of argument as above, we now 
prove that for alH > there exists C{t) finite so that 



Indeed, if we let 
we find that 



_ _ eV^x;'''|U < C{t)e^ (17) 



Yi = -j f 5W^+,ip((l-a)xf + «xf+^'^)tt[y/]cZ5 
Jo Jo 

Jo Jo 

- ^ [DVf3+ei^ - DVfs - e(f^DiVf3]{X^)dt 



Using (16) we find that there exists a finite constant C = C(M) such that 

rt rt 
\\Yf''\\ < C I 

l<i<m 



max < C [ max \\Y^'''\\ds + Ce^ f C{sfds 

l<i<m Jq l<i<m Jq 



and so Gronwall's lemma gives (17) (note here that maxi<j<^ 11^'''' II is finite for all e > so that 

P + €l,eT{c,M)). 

This shows that /3 —>■ X^ is differentiable for all t with first derivative V pX^ . We can continue 
in the same spirit to show that V^xf is differentiable and by induction, we find that /3 —>■ X^ is 
of class C°° in the interior of T(c, M). We next bound uniformly all its derivatives. 

Step 2: P G T(c, M) is analytic for all t > . 

To this end, we can write X^ as a formal series in in a small ball around j3 € T(c, M) 

x?= i: n - ^p)'^^?^ 

(fe)=(fei,...,fe„)eN" l<p<n 

with xf*''"'^^ = X^. Indeed, the coefficients of this series are obtained by differentiating the C°° 
operator X^ and xf^ = {h\ ■ ■ ■ kjy^di^l ■ ■ ■ d^l^X^\,^=p. We write DiVp{X) = E]=iPjD,qj = 

Ylf=i ^jQij where Diqj = YliLi Qi,i+{j-i)D is the decomposition of Diqj as a sum of at most D 
monomials. We denote qij = Yli<,p<di ^i^^ l^j € {1, . . . ,m}. Moreover, we have -qj = V\j/D]- 
Plugging this formal series into 

dX^ = dSt - ^DVr,{X^)dt 
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we find that xf'^i ' • satisfy, for X^fcj > 1, the following induction relation 



dX, 



iD 



=1 y-'*3i fcP-fc _, ,.,^i<P<dj 
A/p=l l^r—l^r J-r=[j/i3] 



(fcP)7^(fc)Vp 



(18) 



z=i 



Above, the sum over indices k such that kr — lr=[i/D] = ~1 is simply empty. Using the convexity 



of V^, we now get a uniform bound by considering X^^Kx^'''^ = YlT=i -^t''^''^ ■-^t''^''' j 



dx['\x^^<-Y^Y. E 



fcr-lr=[i/D] 



l<P<dji l<P<dji 



EE*E E 

j = l 1=1 p=l kP^kr 

(fcp)^(fc)Vp 



dt 



-cx[^\xf'ht 



(19) 



where we have used that by convexity, for any operators X, Z bounded by M, 



Y,Y.^djDiV{X)^Zj).Zi > cZ.Z. 
i=i j=i 

Note that in the right hand side of (19), all the X^^'"''^^^ are such that ^ /^f < ^^j- Hence, we 
can deduce by induction that A^^^ := maxi<j<OT ||oo is bounded for all (k) and uniformly on 

compact sets of the time variable t. Indeed, we proved it in Theorem 2.2 for (k) = (0, . . . , 0). Let 
us assume it is true for (A;) with Y^ki < K — 1 and let us prove it remains true; we simply use 



l<p<dji 



l<p<dij ^<P<dij 



Choosing B such that 



m nD 

2^"'EE(i + i/5*iH-^ E i-c<o 

j=li=l J2k?<kr 
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allows to bound A^''^ uniformly on compact sets by our induction hypothesis. We next show that this 

(k) I (fe) Tk) 

bound can be taken uniformly on the time variable. To this end we first consider = \l Xl .XI 
and deduce from (19) that 

nD 

dYl'^ < E n ^Pdt 

n 

+mD^\Pi\ Yl n APdi-cY^^^dt (20) 

(feP)^(fe)Vp 

where we eventually added terms (by taking dij = D) which can be done if we assume that > 1 
which we can always do. Therefore, for ^ fcj > 1 , since we obtain the bound 



rt 



A^J'^'Us. (21) 



(fcp)^(fe)Vp 

Note that the right hand side of (21) depends only on A^''^ for < — 1 (since < k,. for 

all r but (A;P) 7^ (fc))- Since ^^''^ is uniformly bounded (as we proved that X^ is uniformly bounded 
for /3 G T{c,M) and c > 0), we deduce by induction that A^^^ := sup^>o ^t'^'* is finite and satisfy 
the induction bound, with = Y17=i 1/^*1' 

n 

< mDc-'Y E n ^^"'^ 

i=l v^D _i . l<p<D 

+mDc-^\P\ Y n ^^''"^ (22) 

(fcp)5^(fc)Vp 

We rewrite this inequality, since A^'^^ is obviously finite, as 



-1 



- l + mL>2c-i|/5|(^(o))i3-i\Z^ ii 

"^'^ z=l Y^i? t-P— A- -1 - l<p<-D 

+mDc-^\P\ Y n ^^^^^^l (23) 

,p_, l<p<D 
Z^p=l "^i — — ^ — 

where we added the term mDc~^\P\D{A^^^)^'~^A^''^ to the last sum. We now want to show that 
there is a finite C = C{p) so that 
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To this end we borrow the idea of majorizing sequences as developed by Cartan [9], chapter VII. 
It goes as follows here. We consider the polynomial in one variable given by 

V^{x) = Ckx^ with C :-- '^^'^ ^ 



l + mD2c-i|/5|(^(o))D-i- 
and the the equation 

X-XQ = Vk{x) - V\f3\{xo) 

with xo = A^^\ We claim that for k in a neighborhood of |/3|, the solution in the neighborhood 
of xq of this equation is analytic in k. Indeed, by the implicit function theorem, we only need to 
check that 

which is always true since A^^^ is finite. This implies that there exists a finite constant C = C(|/3|) 
such that 

\d'^x,L=\s\\<k\CK 



But now, x'^'^) := (A;!) ^d^XK\K=\i3\ satisfies the induction relation 

xW = c Yl n ^^''^+cm Yl n ^^"'^ 

Y-D u.-u_Tl<i<D ,._ul<i<D 

which implies 



(1 - C|/3|L>x^-^)xW = C Y n x^'''^ + C\l3\ Y n 

ED u,_u_'\ l<i<D u.—u u.^h.l<i<D 

Since 1 - C\P\Dxq~'^ = 1/(1 + mDc-^\P\DxQ~'^) > 0, we conclude by induction that x^*^) > for 
all k (note that x(°) = A^^^ > 0). But then, comparing (22) and the above inequality, we also prove 
by induction that 

which therefore gives the desired bound for the A^'^K Hence we have proved (15). 
Step 3: Convergence in law of of t ^ {X^''^)(^j^^ for all (fc). 
As we have noticed above, the equations for X^'^'> are of the form 

with some polynomial functions P(j^y Therefore, if we denote X^ the solution of this equation 
starting from X^^ = Z^^^ , we get from the convexity of V and by induction over ^ ki that 

sup ||xf)'^-xf)'°||oo<C(i^)e-* 

with some finite constant C{K). Hence, we can start from Z^^^ = X^^^ to see that since {X^^'^^\^^^n 
has the same law as [X^^'^^i^^^n, the law of [X^'^'^ , ■,^ki< K) converges. □ 
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We now relate the previous result with the absence of phase transition for the generating 
function of colored planar maps. In [20], the following strong version of Schwinger- Dyson equation 
was considered; it requires that for all polynomials P G C(Xi, . . . jX^), 

T{DiVP) = r (g) T{diP), l<i<m. (24) 

It was shown that if V{Xi, . . .,Xm) = W^(gi,/3i)i<i<n = 5 + Td=m+i l^i'H^ ^^^^^ ^^i^ts a 

unique solution M G C(Xi, . . . , Xm)' under the condition that 

\T{Xi,...Xi^)\<R' 

for all k and some finite R, provided the /?i's are small enough. We denote tw, n > this solution. 
Note that this solution was not a priori the law of non-commuting variables, except in the case 
where V is self-adjoint, but just an element of C{Xi, . . . ,Xm)'- In particular, E. Maurel Segala 
and one of the authors always restricted to polynomials in the letters (Xi, . . . ,X^) and did not 
consider their adjoints. 

Moreover, for all monomial P G C(Xi, . . . 

^^(../3,H<,<J^)=>^feAh<.<„(^')= E n ^^^o((g.,fc.),l<^<n,(P,l)) 

fcl,...,A:„eN" l<j<n ' 

with Mo(((/i, fcj), (P, 1)) the number of planar maps with ki stars of type Q'j for 1 < i < n and one 
star of type P (we refer the reader to [20], section 2, for a complete description of the numbers 
Mo((gi,fci),l<i<n,(P,l))). 
We now claim 

Theorem 3.2. (a) The generating function 

mi<i<neC^^M^,„p,),,,,jP):= Yl n ^-^§f^Mi{q,,ki),l<i<n,{P,l)), 

k\,...,kn l<i<n ' 

which is an absolutely convergent series for Y2^=i 1/^*1 •''i^o,^^ enough, extends analytically in the 
interior of the domain where W^g.^p.^^^^^jXi, ... , Xm) = I YliLi + J2i=i PiU (c, M)-convex 
for some c> and M > Mo(c, \\DV{0).DV{0)\\oo) the constant of Theorem 2.2. The extension of 
•^fe,/3i)i<i<n (-^) *^ ^Quo-l to TWi^q. p.)^^.^ ) restriction to C(Xi, . . . ,X^) of the invariant measure 
of Theorem 2.2. 

(h) Assume that V is {c, M)-convex with M > Mq of Theorem 2.2. The invariant distribution fiy 
of Theorem 2.2 not only satisfies (11) but its strong version in the sense that Ty = IJ'v\c{Xi,...,Xm) 
is such that 

Tv{DiVP) = ry (g) Tv{diP), l<i<m (25) 
for all polynomials P G C{Xi, . . . , X^) . 

Hence, the first point of the above theorem shows that the breaking of analyticity (or phase 
transition) of the map enumeration can not take place when t^(qi,/3i)i<i<„ is (c, M)-convex. 

Proof. By [20], if we consider the case where l^(gi,/3j)i<j<„ is self-adjoint, we know that -^{i}i,/3i)i<i<„ 
is the law of self-adjoint operators which are uniformly bounded hy R {R = R{P) going to 2 as /? 
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goes to zero) when the /?'s are smah enough. As a consequence, -^(gi,/3i)i<i<„ niust coincide with 
^i,(ft)i<i<„ = mPi)i<i<u\c{Xu...,x^) where we put /ii,(ft)i<i<„ ■= ^i^fe,^,)i<,<„ ^i^ce MfeA)i<,<„ 
satisfies (11) with potential W(^q. /? (by Theorem 2.2). In particular, ri (/? )i< < must satisfy 

(25). 

We now show that we can remove the assumption that ^(gi,/3i)i<j<„ is self-adjoint. We denote by 
* the involution {zXi^ . . . XjJ* = zXi^ ■ ■ ■ Xi^ so that is self adjoint iff Wf^q.^p.)^^.^^ = 

^te,ft)i<i<„- ^^^^y^ ^{<ii,Pi)i<i<n in the form i^fe,ft)i<,<„,(?*,/3Di<i<n = 2 E"^' + 

^ ^ PiQi + ^ E/^i'Zj*- We denote by jJ-i^jj^p' the invariant measure of Theorem 2.2 corresponding to 
such a potential. The situation V = V* corresponds to f3[ = Pi- In that case we have shown that 

But recall that if an analytic function of two variables x, y is null on A = {x = y, |x| < e} for 
some e > 0, then this function must vanish on its full domain of analyticity since A is totally 
real. Hence, inside the domain of analyticity of [3,(3' '^i.{8i)i<.i<.,-,-(Fi')i<.i<.n (i'^ domain where 
(A)i<i<n, (/3i)i<i<n TiM)i<i<n,mi<i<n^P) is analytic for all P G C(Xi,. . . ,X^)), (26) is always 
true. We can now remove the artificial parameters (/3^)i<i<n to claim that always satisfy 

(25) on its domain of analyticity. We finally can remove the condition that (3i = \ for 1 < i < m as 
follows. We take V (c, M)-convex with M > Mq. Note first that if we consider = f YULi ^1 + ^ 
for some a > 0, then the result still holds since by uniqueness of the solution to Schwinger-Dyson 
equation, we have the scaling property 

with da the dilatation da]\^^{P) = n{P{^, ■ ■ , ^))- As datJF is always (0, M) convex, ^ I] X? + 
dafiV satisfies the above hypotheses and always satisfies (25). Finally, we 

can let a going to zero since is (c, M)-convex for all a > and so a — > /xy^ (P) is analytic and 
thus continuous when a goes to zero. 

As a conclusion, we have seen that ry^ satisfies (26) on the domain of analyticity which contain 
by Lemma 3.1 all the sets r(c, M), c > 0, M > Mo(c). 

On the other hand, we also have that Ti^p must agree with Mq, the generating function of 
maps, for all (3 small enough since ti,^ satisfies (25) and is the law of bounded operators (and it 
was proved that there is at most one such solution, the generating function of maps, in [20]). We 
hence conclude that (3 Ai0{P) extends analytically to the domain of analyticity of ti^^, which 
contains all P such that (l,/3) belongs to T{c,M) for c > and M > Mo(c). 

□ 

4 Connectivity of the support and properties of associated C* and 
von Neumann algebras 

Throughout the rest of the paper, we shall assume that F is a (c, M) convex potential with M > Mq 
(so that the hypothesis of Theorem 2.2 holds). 

In this section, we denote by fi the unique stationary law for the free stochastic differential 

equation (8) with drift DV, satisfying (11), where Sq is a free Brownian motion. Lastly, Z will 
denote some fixed m-tuple of operators, free from Sg : q > 0, having law fj,, and satisfying ||^||cxd < b. 
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The main results of this section concern properties of the C*-algebra generated by the m-tuple 
Z with the prescibed law /xy. We show that this C*-algebra is exact [28], projectionless and that 
the associated von Neumann algebra has the Haagerup approximation property [16] and admits 
and embedding into the ultrapower of the hyperfinite IIi factor. These properties are shared by 
(and in fact, in most cases, derived from those of) the C*-algebra generated by m semicircular 
systems. 

One of the most interesting open problems in operator algebras is a question due to Connes of 
whether any tracial state has finite approximation in the sense that there exists a norm-bounded 
sequence oi N x N matrices {A^ , . . . ,A^) such that for all *-polynomial function P, 

hm lTV(P(Af,...,0)=r(P). 

N—^00 iV 

When m = 1, this question is settled by Birkhoff's theorem, but the question is still open when 
m>2. We prove that the laws fxy have finite-dimensional approximations. 

4.1 Approximation of Z by elements from C*{Sq : q >Q). 

We first show that, with respect to operator norm, C*{Sq : q>0) e-contains a variable with law /x. 

Corollary 4.1. Within the hypothesis and notations of Theorem 2.2, for any e > 0, there exists a 
Brownian motion Sg : s > free from Z , and elements X' G C*{Ss '■ s >Q), X E: C*(Z, Sg '■ s > 0) 
so that X has the given stationary law fi and \\X — X'\\^ < e. 

Proof. Let Xt and X^ be two solutions to (7) with initial data Xq = 0, Xq = Z. By Theorem 
2.2, Xi G C*{Ss '■ s < t) approximates in operator norm the stationary process Xf with marginal 
distribution /x. Since by Lemma 2.1, Xf G C*{Z,Sq : q>0), we may take X' = Xt, X = Xf for 
large enough t. □ 

Theorem 4.2. Let Z be any m-tuple of b-bounded variables with the unique law ji satisfying (11). 
(a) The algebra C*{Z) has no non-trivial projections, (b) The spectrum of any non- commutative 
*-polynomial P in the m-tuple Z is connected (in the case that P{Z) is normal, this means that 
the support of its spectral measure is connected), (c) If P is any polynomial in Z whose value is 
self-adjoint, then the probability measure given by the law of P{Z) has connected support. 

Proof. The C*-algebra A = C*{Ss,s > 0) can be identified with the C*-algebra generated by 
semicircular operators s{f) : f G L^(M;M) where / i— s{f) denotes the free Gaussian functor [4]. 
It is well known that this C*-algebra is isomorphic to the infinite reduced free product 

A={C[-l,l],ii)*---*{C[-l,l],i^) 

where /j, denotes the semicircular measure. The algebra (C[— 1, l],/i) can be unitally embedded in 
a trace-preserving way into the group C*-algebra C*{Z) = C{T) taken with its canonical group 
trace r. Indeed, if it G C*{Z) denotes the group generator, i.e. u = exp{2m9) G C(T), then u-\-u* 
generates a copy of C[— 1, 1], and the restriction of r to C*{u + u*) is the arcsinc law. Hence for 
a suitable continuous function /, f(u + u*) has as its distribution the semicircle law, and we can 
embed (C[— 1, 1],//) by sending its generator, multiplication by x, to f{u + u*). It follows that 

A C C*(Z) *red • • • *red C*(Z) ^ C:^^{¥^). 
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By the results of [24] (see also [17] for a random matrix proof), C*(Foo) has no non-trivial projec- 
tions. Thus A has no non-trivial projections. 

Suppose now that Y £ A. Then the spectrum (t(Y) must be connected. We sketch the argument, 
which can be found in standard C*-algebra literature (see e.g. [6, Proposition 4.6.2 on p. 28]). 
If cr(y) = Ki U K2 with if 1 fl K2 = and both if 1,^2 non-empty, for any contour 7 C C that 
contains Ki but not if 2 and does not intersect cr{Y), the integral 



belongs to A, being a norm limit of Riemann sums. Moreover, = E, E{1 - E) = hut E ^ 0, 
E^\. But -I- - 1 is invertible since (E + E* - Vf = (E - E*){E - E*)* > 1, and thus 



is a nonzero self-adjoint projection (since (1 — E)E = E*{1 — E*) = 0). Clearly, P 7^ and 
(1 — E)P = so P 7^ 1. Hence A would have a non-trivial projection, a contradiction. 

Since G A, it follows that the spectrum of any *-polynomial Y in this m-tuple is connected. 
Since Z and Xf have the same laws at all t, C*{Z) and C*{Xf) are isomorphic. So if the support of 
the spectrum of P{Z) is disconnected, then the support of the spectrum of P{Xf) is disconnected, 
for all t. Because X^ converges to Xf in operator norm as t ^ 00, we find that also spectrum of 
P{X^) must be disconnected for large enough t. Since the trace-state on C*{Ss : s > 0) is faithful, 
it follows that for some t large enough, C*{Ss ■ s > 0) contains an element with disconnected 
spectrum. But as we saw before, this is impossible. Thus we have proved (b). 

Because any non-trivial projection has disconnected spectrum, (a) follows. 

Lastly, as the GNS construction for the trace-state on C*{Z) is faithful (by construction of 
C*{Z)) and the state is tracial, it follows that the trace vector in the representation is cyclic for 
both C*{Z) and its commutant. But this implies that the vector is also separating, so that the 
trace-state on C*{Z) is faithful. Thus the support of the law of any self-adjoint operator X € C*{Z) 
is exactly its spectrum. This implies (c). □ 

4.2 Exactness and the Haagerup property. 

We recall that a C*-algebra A is called exact (cf. [28] and references therein) if there exists a 
faithful *-representation ir : A ^ B{H) with the following property. For any finite subset F C A 
and any e > there exists a finite-dimensional matrix algebra D and unital completely positive 
maps 9 : A^ D, r] : D ^ B(H) so that 



It turns out that this property is equivalent to the statement that the functor (8)min^ of taking the 
minimal tensor product with A is exact (i.e., takes exact sequences to exact sequences). Exactness 
is an important approximation property for a C*-algebra. 

Another important approximation property, this time for a von Neumann algebra, is the 
Haagerup property. A von Neumann algebra M with a trace r is said to have the Haagerup 
property [16] if there exists a sequence of completely-positive maps : M — M, which are unital 




P = E{E + E* -1) 



-1 



Ue{x))-7r{x)\\^<e, 



^xeF. 
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and trace-preserving, so that the associated maps : L^{M,t) — L^{M,t) are compact and 
converge to 1 strongly: 

ll^n(a^) ~ x\\2 — >■ 0, n — > OO 

for all X G M. For a discrete group von Neumann algebra L{T), the Haagerup property is equivalent 
to Gromov's a-T-menability of the group T (i.e., to the existence of a cocyclc c : F — > with values 
in some unitary representation H of F, so that the map 7 1-^ ||c(7)||h is proper). As was shown by 
Haagerup, free groups have the Haagerup property. 

Theorem 4.3. Let Z be any b-bounded m-tuple with the unique law /x satisfying (11). Then C*{Z) 
is exact and W*{Z) has the Haagerup approximation property. 

Proof. As in Corollary 4.1, let Xf be the solution to the free SDE (8) starting with Z and be 
the solution starting with zero. Thus xf e C*{Z,Sq : g G [0, +00)) and Xt G C*{Sg : g G [0, +00)), 
where Sq is a free Brownian motion, free from C*{Z). Moreover, we have by Theorem 2.2 that 

- Xj°||oo as t ^ 00. 

Let E : C*{Z, Sq : q e [0, +00)) — > C*{Sq : q e [0, +00)) be the conditional expectation coming 
from the fact that Z and Sq,q G [0, +00) are freely independent. Then 

\\E{Zt - Xt)\\oo < \\Zt - XtWoo ^ 0, (27) 

so that, since E{Xt) = Xt, 

\\E{Zt) - XtWoo 0. 

Since n is stationary, C*{Xf) ^ C*{Z); let vrt : C*{Z) C*(Xf ) be this isomorphism. Let 
H = L^{W*{Z,Sq : q G [0,+oo))). Then as a module over W*(Z), H is infinite-dimensional, 
and as a module over W*{X^), it is at most infinite-dimensional. Hence there exists an injective 
unital *-homomorphism Gt : B{H) B{H) with the property that @{Trt{x)) = x for alH > and 
X G W*{Z) (and thus for x € C*{Z)). 

Let now F C C*(Z) be a finite subset and e > 0. Then one can find *-polynomials Px : x & F 
so that 

\\Px{Z)-x\\oo <e/9. 
Hence for t sufficiently large, we may assume that 

||i?(P,(Xf))-P,(Xf)||oo<£/9. 

Since -KtiPxiZ)) = Px{Xf), we conclude that 

\\E{Trt{x)) - 7rt(x)||oo < e/3, G F. 

Now for all x & F, E{'Kt{x)) G C*{Sq : q G [0,+oo)), and hence we can find a finite-dimensional 
C*-algebra A and unital completely positive maps rj : C*{Sq : g G [0, -|-oo)) ^ A, ip : A ^ B{H), 
so that 

IIV" o »7(y) - J/lloo < e/3 
for all y G 7rt(F). Consider now the unital completely positive maps 

a = r}oEoTrt:C*{Z)^ A 
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and 

P = eotp : A ^ B{H). 

Then 

\\l3oa{x) -x\\oo = Wrjiirtix)) -7rt(x)||oo < £ 

for all X G F. Thus C*{Z) is exact. 

We now turn to the Haagerup property, where we adapt a proof from [7]. Consider the map 
$t : W*{Xf) W*{Xf) which is obtained as the composition 

= Ew*(xf) ° *t ° Ew*{Sg:qe[Q,+oo)), 

where are unital trace-preserving completely positive maps on W*{Sq : q G [0, +00)) so that 
are compact on and — x\\2 — >■ for all x G W*{Sg : q G [0,+cx))). Then are 

unital trace-preserving completely-positive maps on W*{Xf ) = W*{Z), and because of (27), one 
has that ||$t(a;) — x\\2 — > for all x G W*{Z). On the other hand, each <&t, viewed as a map on 
L?{W*{Z,Sq : q G [0, -|-oo))), is compact (since is compact), and therefore the restriction of $t 
to L'^{W*{Xf)) is also compact. Thus W*{Z) has the Haagerup property. □ 

4.3 Finite dimensional approximation. 

4.3.1 embeddability for self-adjoint potentials. 

In this section, we improve on the results of [20] by showing that if we set 

with 

Aiv(P) = ^TV(P(Ai,...,A^)) 

and 

then p.y converges towards ry for any self-adjoint locally strictly convex potential V. In [20], a 
similar result was proved when V = cX.X + W with W a 'small' enough polynomial. 

Theorem 4.4. For all c > 0, there exists Bq < 00 and Mq < 00 so that for any self-adjoint 
polynomial V which is {c, M)-convex with M > Mq there exists a unique law t of m self-adjoint 
variables such that for all i E {!,..., m}, all polynomial P, 

T ® T{diP) = T{DiVP) (28) 

and such that T{Xf^) < B^^. Moreover, Jly converges towards r and therefore r has finite approx- 
imation. In particular, if Z has law t, then W*{Z) can be embedded into the ultrapower of the 
hyperfinite IIi factor. 
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Proof. We can follow the lines of [20] Theorem 3.5 to see that d/x^ (Ai, . . . , Ajn)) has a density 
with respect to the Gaussian law Yl e dXi/Zf^ which is log-concave. This insures that we can 
use the Brascamp Lieb inequality which in turn allows us to show that the random matrices under 
the above Gibbs measures stay bounded in norm by some Bq <^ M with overwhelming probability. 
As a consequence, we can perform an infinitesimal change of variables Xi ^ Xi + N~^Pi with 
Pi a self-adjoint polynomial in ||A||oo < -Bo and the null function outside a ball of radius strictly 
smaller than M. This shows that any almost sure limit points fij^ under are laws of variables 
bounded by Bq and satisfy (28) with the potential V (see [20] for details). Now, if they satisfy (28) 
they also satisfy (11) (take P = DiP and sum the equalities) and so by Theorem 2.2, there exists 
at most one such solution fiy = tv- Thus, /tjv converges almost surely and therefore in expectation 
towards ry. Hence ry has finite approximation. □ 

4.3.2 R'^ embeddability for non-self adjoint potentials. 

One can give a proof of embeddability of W*{Z) into the ultrapower of the hyperfinite IIi factor, 
based directly on Corollary 4.1. This proof works for arbitrary (c, M)-convex polynomials (without 
the self-adjoint assumption). Indeed, because of this corollary and with its notations, if e > is 
given, and F is a finite collection of *-polynomials, then there exists an X' G C*{Sq : q > 0) with 
the property that 

\T{P{Z))-T{P{X'))\<e, VPGF 

(here r denotes the free product trace-state on C*(Z,Sq : q > 0)). Now, C*(Sq : g > 0) is 
generated by an infinite free semicircular family Sj : j = 1,2,.... One can clearly assume that 
in fact X' G C* (5*1, ... , Sk) for some large enough K (since one can replace X' with Ek{X') for 
K large enough, where Ek '■ C*{Si : i = 1,2...) C*{Si, . . . , Sk) is the canonical conditional 
expectation). Thus, by approximating X' with a polynomial in Si, ... , Sk, we may assume that 
X'l^ is a polynomial Qk in ^i, . . . , Sk- Since Si, ... , Sk are free semicircular variables, their law 
has finite approximations. Thus for N sufficiently large, one can find a iiT-tuple oi N x N self- 
adjoint matrices ^41, . . . , Ak whose law approximates that of 5*1, ... , Sk so well that the m-tuple 
of matrices B = {Bi, . . . , Bm) = (Qi(^i, • • • , Ak), . . . , Qm{Ai, . . . , Ak)) would have the property 
that 

|r(P(Z)) - 1tV(P(B))| < VPGP. 

5 Free entropy 

For the remainder of the paper we shall assume that V is (c, M)-convex and self-adjoint. 

In this section we show that for tracial state with conjugate variables given as the cyclic gradient 
of a self-adjoint (c, M)-convex potential the microstate entropy is the same whether it is defined 
by a limsup or a liminf. 

Theorem 5.1. Let c > and V be a self-adjoint (c, M)-convex potential with M > Mq. Let t = Ty 
be as in Theorem 4-4- -^ei T{T,e,R,k) be the microstates 

r(r, e,R,k) = {Xi,--- ,Xm: \\Xi\\^ < R,\j^Tr{P{Xi, ■■■ , X^)) - r(P)| < e, 

for all monomial of degree less than k} 
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and let vol denote the volume on the space of m N x N Hermitian matrices. 
Then 

1 TYl 

x(r) = lim sup lim sup —k log uoZ(r(r, e, i?, k)) + — log N 

eiO,k,R-^<x> N-^oo 2 

1 Tn 
= lim inf lim inf — ^ log vol{r{T, e, R, k)) + — log N 

eiO,k,R-*oo W— »oo N 2 

and x(r) > — oo. 
Proof. Note that 

V0l{r{T,€,R,k)) = J lAi,-,Arr.&T(r,e,R,k)d^l---dAm (29) 

ir(r,e,R,fc) 

V{)^(r(T, e, R, k)) J e-^'^''(^)d^i ■■■dAm (30) 



with yu]^^ the Gibbs measure considered in the proof of Lemma 4.4. Here, we used the notation 
A{N,e) « B(N,e) when 

1 , ^(A^,e) 
lim lim —77 log — r = 1. 
elo JVTO ^ B(iV, e) 

Since iiy{T{T,e,R,k)) ^ 1 by the proof of Theorem 2.2, we only need to estimate the quantity 
/ e~^'''^'^(^)o?yli • • • dAm- To do that we write V = W + ^X.X with W a (0, M)-convex potential 
so that 

a,i,log/e-^^'^"(*'^+f^-^)d^,...d^^ = -A,V|x.x(^) 

By the proof of Theorem 4.4, we also see that for all t G [0, 1] 

lim Jlm+^x.xi^) = nw+^x.x{W) 

(note that tW + ^X.X stays self-adjoint (c, M)-convex). Since everything stays bounded, and since 
the limit 

hm log / e-^^^'^(f^-^)d^i ... dA^ + log N 

W— >oo 1\ J Z 

is a finite constant -F(c), we conclude by bounded convergence theorem that 

Jirn^ ^ log I e-^''^"(^+t^-^)d^i . . . + ^ logiV = F(c) - j^' w+|x.x(VF)dt 
which shows by (30) that % can be defined either by limsup or liminf and both are equal to 

x(t) = F(c) - r^^^ix.xmdt + t(I^) 

Corollary 5.2. Let r 6e as m Theorem 4-4- Then the von Neumann algebra generated by an m- 
tuple Z with law t is a factor, does not have property T, is prime and has no Cartan subalgerbas. 

Here we of course use the fact that xi^) > ^^'^ that for an m-tuplc Z, xi^) > implies 
that W*{Z) is a factor, non-F, has no Cartan subalgebras [26] and is prime [15]. 
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6 Norm convergence 



In section 4.3 it was shown that solutions to the Schwinger-Dyson equations (28) are weak hmits 
of finite dimensional approximations. Namely if we let 



mi 



with 

iiN{P) = ^TT{P{Au...,Am)) 

and 

dtx^{A^,...,Am)) = ^l||^,||<Me-^'^(^(^^'-'^'"»dAi...dA 

we saw that ]Ji^{P) converges to the unique solution r of 

Tv{DiVP)=Tv®Tv{diP) 

for alH G {1, . . . , m} and P G C(Xi, . . . , X^). We here show that this convergence holds in norm 
but for simplicity restrict ourselves to potentials which are uniformly convex (and not only locally 
convex) . 

Lemma 6.1. Let (Xi,...,X^) he a m-tuple of non- commutative variables with law ry and let 

[X^ , . . . , X.I^J be random matrices with law Hy . Assume that V is self-adjoint {c,oo)- convex for 
some c > 0. Then, for any polynomial P G C{Xi, . . . , X^), 

hm ||P(Xf , . . . ,X^)||oo = \\P{Xi, . . .,Xm)\\oo a.s. 

This result generalizes the work of Haagerup and S. Thorjornsen [14] where it was proved for 
V{Xi, . . . , Xm) = 5 YliLi Xf- Our result actually relies on theirs. 

Proof. The idea is to use the approximation by processes, and hence the fact that processes are 
well approximated by polynomials of independent Wigner matrices and then use [14] to conclude 
that the norm of the latter converge to the limit. 

Step 1: Matrix valued diffusions and convergence to the stationary process. 

Consider the diffusion with values in the set of Hermitian matrices 

dXf = dijf - ^-DV{X^^^)dt (31) 

with X^'^ = Z. Here {H^ , t > 0) is a m-dimensional Hermitian Brownian motion. In other words, 
i?/^ are a set of m independent matrix valued process whose matrix entries are given by 

rNn„ i\ Bkl{t) + iBki{t) fjNtu i\ uN/u u\ ^kk{i) 



HI' (k, I) = k < I, HI' (l, k) = HI' (k, I), HI' {k, k) 



2N VN 

where the B, B are independent standard Brownian motions. 

A strong solution to (31) exists up to a possible time of explosion (when DV would stop being 

Lipschitz eventually) since this equation can be seen as a system of classical stochastic differential 
equation with N{N + l)/2 equations driven by independent Brownian motions with a polynomial 
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drift. Since this drift derives from a strictly convex potential, it is well known that the time of 
explosion is almost surely infinite (which can also be show by the arguments of the proof of Lemma 
2.1) 

Now, let us consider two solutions X^'^ and X^'^ starting from Z and Z' respectively. Then, 
we have 

- '^') = -\{DV{X^^^) - DV{X^'^'))dt. 
Hence, we can apply exactly the same arguments than in the proof of Theorem 2.2 to conclude that 

ll^t ~ lloo ^ e \\Zj — Zj Woo- 

If we take Zn to be random with law jJLy we get a stationary process so that 

||^jv,z _ x^'^W^ < e~'^*||Zjv||oo- (32) 

Now, according to Brascamp-Lieb inequality (see its application on our particular case in [21]) 

Mydl^Tvlloo >x) <e-«('^)^(---°W) 

with some xo{c) < oo and ao(c) > 0. Therefore, (32) shows that for all t>0, 

P(||xf '°||oo > (e-^* + 1)K) < 2e-'^('^)^(-^-^o(^)) (33) 

since both Zpf and X^'^ have law ^y. 

Step 2: Uniform hounds on X^'^. In this section we want to show that we can also control 
uniformly the norm of X^^ = X^'^ for t in a compact set. To do that let us remind that 

dX^.X^ = 2X^ AH^ - DV{X^).X^dt + 2dt 

< 2X^.dH^ - cX^.X^dt + 2dt - 2DV{Qi).X^dt 

< 2X^.dH^ - ^cX^.X^dt + Cdt 

with C = 2+ "^DV {{)) .DV {{)) . Therefore, for any p > 0, since X^ .X^ is a non negative matrix for 
all s > 0, Ito's calculus yield 

T,{{X^.X^Y) < 2p j\r{{X^ .X^Y-\X^ .dH^' -\cX^ .X^ dt + Cdt)) 
+2pN-' r Tr((Xf .Xf )^)Tr((Xf .Xf )P-*^-i)dt 

+2pN-^ r Tr((Xf .X,^)^X,^)Tr((Xf .Xf )^'-'^-2xf )dt (34) 
Now, by Burkholder-Davis-Gundy inequalities and Chebychev's inequality, there is a finite constant 
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A such that for all e > 0, 
Jo 



sup 

, s<T 



> eN2 < 



A 



iiV2 



E 



1 riv((xf 





dt 



< 



AT 



where we used (33) and chose K = max{2xo(c), 1}, p < a{c)NK. Therefore, if we choose e so that 
e = {K+ VfPT, we see that if we set 

A{N,T)= fl (sup f Tr{{X^ .X^Y-\X^ .dH^) <{K + lfPTm\ 
p<a(c)NK J 



p<a{c)NK 

then by Borel Cantelli's Lemma, 



limsupP(^(iV,r)) = 1. 

7V,T^oo 



Let now restrict ourselves to the set r\T>To f^N>No ^{^i T). We let 

Bip,N,s) := N-'Tr{{X^.X^r) 
and observe that these non negative real numbers obey the relation 

B{p,N,s)<B{q,N,sf^, q > p. 
We first control B{po, N, s) by (34) which yields for s > 0, 



o 1 f PQ-J- 

B{po, N, s) < 2po{K + 1)2P0 Ar-2 s V Tq + / {-2cpB{po, N, t) + 2po(2po - l)B{po, N, t) po )dt 

Jo 

<2po(i^+l)'^"iV-2sVTo+ / {~2cpB{po,N,t) + 2po{2po-l)B{po,N,t) ™ )dt 

Jo 

< 2po{K + l)2foiV-^s V To + { [ - 2cpo + a2poi2po - 1)] B{po, N, t) 



+ 



/ po - l \Po-i _ /po^\Po 
V poa J V vna J i 



poa 



[dt 



where wc used that x «) < ax + ( ^°^ )Po-i _ a( ^^^ )^° for all a > and po ^ 1- Choosing a so 
that a(2|)o + 4po(]3o — 1)) = c^^o we thus have found a finite constant C{po, K) such that 

B{pQ, N, s) < C{po, K)s V To - cpo I B{po, N, t)dt 

Jo 

which shows that B{pq,N,s) < m.ax{C{po, K)/cpo,C{po, K)Tq} is uniformly bounded. 
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We now bound B{p,N,s) for p > po and to this end replace in (34) all B{q,N,s),q < po by 

PQ 

B{po, N,s) 1 . We show by induction over p>po that B{p, N, s) < CpCg with Co a finite constant 
depending on po and Tq, for all s > and p < a{c)KN. Here Cp denotes the Catalan numbers. 
Indeed, this is satisfied for p < po and then (34) implies that 

B{p,N,s) <2p{K + l)^PN-hvTo+ [ {-2cpB{p, N,t) + 2pB{p, N,t)'^)dt 

Jo 

P-l i-s 



dt 



where we used that \^Tr{{X.X)''X)\ < ^Ti{{X.X)''+^). Because 1 < Q < 4Cfe_i and since 
Ylk=o ^kCp-i-k = Cp for all k and p we conclude that 

B{p,N,s) < 2p{K + lfPN-^s\/To+ [ {-cpB{p, N,t) + 2pB{p, N,t)^)dt 

Jo 

+WpCP-^CpS. 

< 2p{K + ifPN-^s V ^0 + (^]^cpB{j),N,t) + 0^ ^dt + lOpC^'^CpS. 

Thus, if N is large enough so that N~'2Tq < 1, we get that 

B{p,N,s) < (^{K + lfPN-'2 + (^^J + 10Cl%^ 

with /(To) = To if s < Tq and /(To) = | if s > Tq. Note here that we used the fact that we have a 
negative drift growing linearly with p to cancel the multiplication by p. We finally choose Cq large 
enough so that 

M^(^(K + l)^PN--2 + (^^J + 10CP-%^<CP 

which we can always do. 

Hence we have proved that 

sup .Xf lloo < min sup{NB{p, N, s)f/P < 2C'{po) < oo a.s. 

t>0 p<a{c)KN t>0 

In other words, we have proved that HX^^Hoo is uniformly bounded almost surely. 

Step 3: Convergence of the norm of P{X^'^) as N goes to infinity. To this end remark that 
since Xt is always uniformly bounded by M we can always assume V is uniformly bounded 
and with uniformly Lipschitz cyclic gradient (this amounts to change V outside a place that the 
diffusion does not see). We let V be equal to V on operators with norm bounded by M and have 
uniformly Lipschitz cyclic gradient. For instance, we take V{Xi, ... , X^) = V{f{Xi), ... , f{Xm)) 
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with f{x) = X on \x\ < M, f{x) = x{l + (|x| — M)^)~^ if |x| > M (since the later is twice 
continuously difFerentiable with uniformly bounded derivatives). 
Now, by definition if we let 



riX,S)t = St-^J^ DV{X,)ds, 

then Xt can be expressed as an iterate 

Xt = (l>M{X,S)t = (l)M{.,Snx)t 

for all integer numbers n and M greater than the uniform norm on Xt,t > 0. On the other hand, 
for two operator valued processes {X, Y) 



UM{X,S)t-cPM{Y,S)t\\oc<\\\DV\\cJ^ 



and so we get that 

\\xt-cp{.,sns)t\\o. = mx,s)t-cpir-\.,s){s),s)t\ 



< \\DV\\c f\\Xs-r-HS,.){SUoods 
Jo 



< C- — z — - — TT-sup 
2"(n-l)! „<^" 



LDF Kir „ ~ 
< C " , „ DV 
2"(n-l)!" 



oo 



for all n G N. 

We next want to show that the norm of P(X^'^) converges with overwhelming probability for 
any P G C{Xi, . . . , X^). To do that we approximate X^''^ by (pMi-, -f^^)"(0)t with as above 
We claim that 

lim lim lim sup ||0m(-, H^TiO)t - ||oo = a.s. 

AT— *oo M-+0O n— »oo j< 

Indeed, 4>m{-, H^) is a contraction for all M finite and {X^ , t < T) is its unique fixed point as long 
as {Xf^ ,t < T) stays uniformly bounded by M. Since we have seen that almost surely {Xf^ ,t < T) 
is uniformly bounded, the statement follows. 

Step 4: Convergence of the norm of(j)M{-, H^). For all M, 4>m{-, H'^)"'{0)t can be approximated 
uniformly by a polynomial function of on sup^^^ ||i/^||oo < L} for some L finite, which happens 
with probability one for some sufficiently large L. We can thus use [14] to conclude that the norm 
of any polynomial in (f)M{-,H^)^{Q)t converges to its analog with replaced by S. 

Step 5: Conclusion. We have proved that for all n and M 

TN\ni 



hm ||P(0M(.,//^'r(O)t)||oo = ||P(</>M(.,5r(0)t)||oo a.s. 

iv— »oo 



Since X^ is uniformly bounded, for N large enough. 



lim lim sup ||(/.m(., H'^T{^)t - ||oo = a.s. 
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implies 

lim lim snp\\P{4>M{.,H^T{0)t)-P{Xt^)\\oo = a.s. 

M— >oo n^oo ^^rp 

And finally we have with overwhelming probability (where with stationary law such that 
||-^Ar|| < K) for N large enough 

Let e > be fixed. We fix t so that e~'^'^K = e. With overwhelming probability, since X^'^ has the 
same law than Z, \\X^'^\\^ < K and so also ||X/^||oo < e + K. Hence, we have for any polynomial, 

hm - )|U < C{K)e a.s. 

iV— >oo 

Now, we take M, n large enough (greater than some finite random integers) to assure that for N 
large enough so that 

\\P{<f>M{., H''nO)t) - P(Xf )||oo < 6/3 

and finally 

\\\P{4>M{.,H''r{0)t)\\oo - ||P(.^M(.,^f (0)t)||oo| < e/3. 
We have already seen that 

\\P{cPM{.,SnO)t) - P{Xt)\\oo < e/3 

Thus, we have proved 

|||p(xf'^)IU-|imf)|U|<c'(K)6. 

This completes the proof since X^'^ (resp. Xf ) has the same law that Zn (resp. Z). □ 
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